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We bosonize the long-wavelength excitations of interacting 
fermions in arbitrary dimension by directly applying a suit- 
able Hubbard- Stratonowich transformation to the Grassman- 
nian generating functional of the fermionic correlation func- 
tions. With this technique we derive a surprisingly simple ex- 
pression for the single-particle Greens- function, which is valid 
for arbitrary interaction strength and can describe Fermi- as 
well as Luttinger liquids. Our approach sheds further light 
on the relation between bosonization and the random-phase 
approximation, and enables us to study screening in a non- 
perturbative way. 
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Several groups jl]-[I[ have recently constructed 
bosonization rules for interacting fermions in dimensions 
d > 1. The strategy adopted in these works follows 
closely the usual bosonization of one-dimensional systems 
H , and is based on the observation that suitably defined 
local density operators approximately obey bosonic com- 
mutation relations in the Hilbert space of states with 
wave- vectors close to the Fermi surface. An alternative 
method to bosonize interacting fermions without calcu- 
lating commutators is based on functional integration 
HQ . In this letter we shall further develop this approach, 
and show that it is in many respects more powerful than 
the usual operator bosonization. 

Consider a system of interacting fermions on a d- 
dimensional hypercube with volume V. The hamiltonian 
is given by H = H Q + H mt , with H = £) k e k c k c k and 

Hi nt = ^jEqkk' /q k ' a k+ q 6 k'-q 5 k'Ck, where c k annihi- 
lates an electron with wave- vector k, and is some ar- 
bitrary energy dispersion ||. We assume that the de- 
grees of freedom far away from the Fermi surface have 
been integrated out, so that k and k' are restricted to 
a shell around the Fermi surface, with radial thickness 
small compared with the Fermi wave- vector kp. Thus, 



E k © Q (k)c k Ck+q, where 9 Q (k) = 1 if k e K a , 
and zero otherwise. Assuming that the variations of / kk 
are negligible if k and k' are restricted to given patches, 
we may introduce the coarse-grained interaction function 
yW =< e Q (k)6 Q '(k')/q k ' > where < . . . > denotes 
averaging with respect to k and k'. 

To bosonize H, we consider the imaginary-time corre- 
lation function 

= y f dre-^ < i%(r)iK q (0) > , (1) 

where /3 = l/T is the inverse temperature and q = 
[q, iui m ] denotes wave- vector q and bosonic Matsubara 
frequency io m — ^TtmT. Below it will become evident 
that bosonization of H is equivalent to the calculation 
of a functional S e f / {p} of a complex field p q , such that 
Eq.[l] can be written as a bosonic functional integral 

, p jV{p\p a q p't q e- s 'ff^ 



/ kk and £k differ from the bare jparameters by finite 



renormalizations. Following Refs. [1-|4[], we partition the 
Fermi surface into disjoint patches of volume A d_1 . The 
precise shape of the patches is irrelevant as long as we 
are interested in fluctuations with wave- vectors |q| <C A, 
and the interaction / kk is dominated by momentum 
transfers |q| -C kp. Throughout this work we shall as- 
sume that these restrictions are satisfied. Let us intro- 
duce a label a to enumerate the patches in some arbi- 
trary ordering, and denote by K a the set of k-points in 
patch a. For each a we define local density operators 
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V 



f{p} 



(2) 



JV{p}e- 

We now derive S e ff {p} in arbitrary dimension. Starting 
point is the representation of H q a as a functional integral 
over Grassmann fields c k § 

e -Si{c}-Si{c} 
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(3jV{c}R«R% 



V 



-S 1 {c}-S 4 {c} 



(3) 



with Si{c\ = P^2 k [-iu n + €k]c' k Ck, and S 4 {c} = 
|E 3 Ea a ^-,/q M '<' Here Q n = n(2n + 1)T are 
fermionic Matsubara frequencies, £k = £k — P is the en- 
ergy measured relative to the chemical potential p of the 
interacting system, and R q = J2k ® a (kO c ! c fc+<?- We now 
decouple S4 {c} by means of a Hubbard-Stratonowich 
transformation that involves two auxiliary fields tfr™ and 
Pq ©■ Completing squares and using the invariance of 
the integration measure with respect to shifts ||, it is 
not difficult to show that 



S4,{c} _ 



V{p}T> {</>}< 



-S 2 {c,4>}-S 3 {4>,p}-S i7lt {p} 



(4) 



s 2 {c,4>) =J2 c I v m' c « > ^=*X> a ( k )$ 1 -fe' (5) 

kk' a 
q a 

^w4eea« . /r'=^/r'- (7) 
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The field p^ is the bosonized local density field, while cj)^ where U^(— qa, qa') = ^IIq q (q), with 
is the field dual to p^. Because the composite fermionic 
field Rg couples to the dual field 0™, the latter must be 
introduced at least at an intermediate step in order to be 
able to integrate the fermions out at the very beginning. 

Because the effective fermionic action is quadratic, we 
may now perform the Grassmann integration, and obtain 



n r'(9) = ^D "( k ) ea '( k + q) 



/(&) ~ /(& +q ) 



(12) 



TTQQ 



(8) 



P J V {p} V {(/)} p«p« qe -s^M-s 3 {4>.p}-s mt {p} 
V J X> {p} X> {</>} e-^™ W-s 3 {0, P }-s lnt {p} 

where Skin {</>} can be expanded as Skin {<f>} — 
EZiSHm, with = ^Tr[G F]". The 

trace is over wave- vector and frequency space, and 
Go and V are infinite matrices, with [Go] kit' = 
<5k,k'^n,n'G (fc), and G (k) = [iw n - £k] _1 - The matrix- 
elements of V are defined in Eq.||. The subscript on 
Skin indicates that this quantity is closely related to 
the bosonized kinetic energy. Performing the trace in 

Sj&{4>] Yields 



qi...q n oi...a„ 

x C/ (n) (giai 



+q„,0Omi+...+m n ,Q 



(9) 



f/W(^...^) = (j) n ^EC 



(A:) 



x Go " 2 (fc + q Pl ) ■ • ■ Gq P « (fc + q Pl + . . . + qp n _ t ) , (10) 

where G%(k) = 9 Q (k)G (fc). We have used the in- 
variance of S^ {</>} under relabeling of the fields to 
symmetrize the vertices with respect to the inter- 
change of any two labels. The sum ^ p is over the 
n\ permutations of n integers, and Pi denotes the im- 
age of i under the permutation. Note that the vertices 
are uniquely determined by £k- In d > 1 or for 
one-dimensional models with non-linear energy disper- 
sion, ((?!«! . . . q n a n ) are rather complicated func- 
tions of all external momenta and frequencies. However, 
in the Tomonaga-Luttinger model (TLM) (i.e. for one- 
dimensional fermions with linearized energy dispersion) 
all with n > 3 vanish identically. This is a direct 
consequence of the closed loop theorem, which is discussed 
and proved in Ref . |l(J . Note that it is necessary to sym- 
metrize the vertices in order to apply this theorem. 
The first term in the expansion of Skin {4 1 } is 

S { k ]im = *E a N a ^, where N° = £ k e«(k)/(£ k ) 
is the number of occupied states in patch a. Here 
/(e) = [e@ e + 1] _1 is the Fermi function. For any fi- 
nite q the contribution from <S^„ cancels in Eq.|| The 
second-order term is 



(ii) 



The higher order vertices describe non-local interactions 
between the fields. In the absence of nesting wave- vectors 
and Van Hove singularities, they are non-singular and 
have a finite limit if all external momenta and frequencies 
are set equal to zero jll). In this limit is diagonal 
in all patch-labels, and is proportional to the (n — 2) nd 
derivative of the density of states at the Fermi surface 
with respect to the chemical potential jlj]. Note that 
in the TLM the density of states is approximated by a 
constant, so that it is also from this point of view clear 
that U^(qi — 0) vanishes for n > 3. However, the closed 
loop theorem is stronger, and guarantees that the higher- 
order vertices vanish even at finite values of the qi . 

At this point we integrate over the field </>. The effective 
action S e ff {p} in Eq.0 is given by 



Seff {p} = S lnt {p} - In 



V{4>}e 



-S kin {<j>}-S 3 {<t>,p} 



(13) 



Performing the integration in Eq|l3| perturbatively, 
we obtain an expansion S e ff{p} = J2n°=i ^ef\ {p}' 
where S^ {p} is of the same form as Eqj9[ but with 
the ^-fields replaced by the p-fields, and the vertices 
ai...q n a n ) replaced by a new set of vertices 
T^ n > {q\a\ . . . q n a n ) . In general contains contribu- 
tions from all [/( m ) with m > n. We now truncate the 
expansion of Skin {0} a t the second order. In the TLM 
this approximation is exact, while in d > 1 it is expected 
to be accurate at high densities and for models where the 
higher-order derivatives of the density of states are small. 
The integration in Eq.O is then Gaussian, and we obtain 
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-qa, qa 



(14) 

where T^ 2 ) is the matrix inverse to in the space 

spanned by the patch indices. Since we are interested 
in long-wavelength fluctuations, we may approximate 
r'- 2 -' by its leading term for small |q|. At high densi- 
ties U^ 2 \— qa, qa 1 ) cx 8 aa , because the sum in Eq.|l2| is 
dominated by momenta k of the order of fcp, so that we 
may set 9 Q (k)e a '(k + q) w 5 aa '@ a (k) in Eq£§ In this 
approximation we obtain 



r(a)(. 



■qa, qa 



(15) 



where v a = V~ 1 dN a /dp, is the local density of states in 
patch a, and f£ = v° • q, with v Q = Vk£k|kei<" Q - Insert- 
ing Eq.||| into Eqjl|, we see that the term proportional 



2 



to iuj m defines the dynamics of the p-field. We now recall 
that in the functional integral for canonically quantized 
bosons the coefficient of the term proportional to —iu> m 
should be precisely [3. Thus, to write our effective action 
in terms of a canonical boson field b q , we should rescale 
the p-field accordingly. This is achieved by substituting 
Pi = (V*/^|)*[e(^+e(-^)&l a 8 ] inEq.0. Ourfi- 
nal result for the bosonized action {b} = S^X {p(b)} 



of the energy dispersion Q a can be calculated exactly 
for a given realization of <j), because the differential equa- 
tion for Q a is first order and can be easily solved 0. 
The physical Greens-function is then obtained by aver- 
aging Q a over the distribution of the field (f>. Within 
the non-interacting boson approximation, this average 
involves a trivial Gaussian integration. The final re- 
sult for the real-space imaginary-time Greens-function is 
G Q (r,r) = G£(r,T)e3 Q ( r - r \ with 



q a. 

+ f3[H kin {b} + H int {b}} , 

q a 
q olol' 

x e(q') + fz' a bfb«) 

+ e(-q') (ffbfb a :i + ]^ a b a : q b a q 



(16) 



(17) 



(18) 



where f® a = \/v a v a f^ a are dimensionless couplings. 
The functional integral for the 6-field is now formally 
identical with a standard bosonic functional integral. 
The corresponding second quantized bosonic hamiltonian 
is therefore H b = H b kin + Hf nt , where H\ m and H b ?t 
are simply obtained by replacing the bosonic fields b q in 

Eqs.|T7|,|l8]by operators b™ satisfying [b^,b^] = 5 aa 'Sq^. 

The resulting H b agrees with the bosonized Hamiltonian 
derived in Refs. |^|| by means of an operator approach. 
From our derivation it is obvious that bosonization is a 
non-perturbative but approximate method in d > 1. The 
non-interacting boson-approximation is only accurate if 
the higher order terms Sj^l, n > 3, are irrelevant. 

H b and S^tt {p} (EqJIJ) contain the same physical 



information. To calculate IT?" , it is convenient to 
parametrize the functional integration in terms of the 
original density field p q . Inserting Eq|l4| into Eq.^ we ob- 
tain IT;"' = [[U^ 1 (q) + fc i r 1 } aa ' where n and / should 
be understood as matrices in the patch labels. This 
is nothing but the random-phase approximation (RPA), 
which is known to be exact in the TLM |l2| . The phys- 
ical density-density correlation function is J2 aa , II" a . 
In d > 1 the summation over the patches gives rise to 
Landau-damping. 

Let us now focus on the single-particle Greens-function 
G a (r, r) associated with patch a |Q . It is well-known 
jnj that the Hubbard-Stratonowich transformation re- 
duces the calculation of G a (r, r) to the problem of cal- 
culating the Greens-function Q a {r, r', r, r'; {</>}) of an ef- 
fective non-interacting system in a dynamical random- 
field <b. The crucial observation is that after linearization 



Q Q (r,r) = R a -S a (r,r) , R a = lim S a (r,r) , (19) 

r,r—*Q 

SM -^v\ — — ■ () 



Here f q RPA — [f^ 1 + Ho (?)]" 1 is again a matrix in the 
patch indices. Eqsjl^ and [2^ can be used to determine 
by direct calculation of the Greens-function if the sys- 
tem is a Fermi- or Luttinger liquid, orperhaps belongs 
even to a different category. In Refs. |pyi4f the Greens- 
function is also found to be of the form G a = G^e® . 
Note, however, that Rcf. [Q and Ref. [0 give differ- 
ent expressions for Q a . Our Q a resembles more the re- 
sult derived by Castellani et al. via Ward identities, al- 
though it is not precisely the same, because their Q a de- 
pends explicitly on an ultraviolet cutoff. In our approach 
all cutoff-dependence is contained in the local density of 
states v a . Below we show that for a quite general class 
of interactions Eq.|2(] depends only on the global den- 
sity of states v = '^2 a v a , and is therefore manifestly 
cutoff-independent. We now discuss Eqs|n] and ^ for 
V,f3 — > oo in some physically interesting cases. A de- 
tailed derivation will be given in a longer publication. 

(ofLuttinger liquid in d = 1. In one dimension Eqs.|l9| 
and |20| correctly reproduce the known results for the 
Greens-function of the TLM. Note that R a and S a {r,r) 
are both logarithmically divergent if /°° goes to a finite 
limit for q — > 0, but Q a (r,r) is finite. It seems that 
the possibility of expressing the interaction dependence 
of Q a (r,r) exclusively in terms of the RPA- interaction 
has not been noticed in the literature on the TLM. 

(b) Quasi-particle residue. If the integral defining R a 
exists, the system is by the usual definition a Fermi 
liquid. The quasi-particle residue Z a for wave- vectors 
k € K a is then Z a = e Ra . From Eq.|o| it is not difficult 
to show that R a is real and negative, so that < Z a < 1. 

(c) Singular interactions. For interactions of the form 
/q Q = g 2 |q| _r 'e _ ' q '/ 9c , r\ > 0, we find by simple power 
counting that R a exists only for r) < 2(d — 1), in agree- 
ment with Ref. ]l5| . If in addition < rj, and if the 
screening wave-vector k = {vg 2 ) 1 ^ 11 is small compared 
with hp, then the ultraviolet cutoff q c is not necessary 
and may be set to infinity, because re acts as a natural 
short-distance cutoff for all wave-vector integrals. It is 
important to stress that we are not a priori assuming 
that the bare interaction is screened. A straight-forward 
calculation gives in this case for the quasi-particle residue 



3 



exp 



r(d, rf) 



2(d-l) -17 



{k/Uf) 



d-l 



(21) 



where the positive numerical constant r(d, rf) can be ex- 
plicitly written down as a one-dimensional integral, and 
remains finite at r\ = 2(d — 1). 

(d) Coulomb interaction in d > 1. This is a special 
case of (c): g 2 = s^e 2 , r\ = d — 1, and q c = 00. Here 
Sd is a numerical constant (S2 = 2tt, S3 = An), k = 
(vsde 2 )' 3 ^ 1 can now be identified with the usual Thomas- 
Fermi screening wave- vector. Note that k <C kp in the 
high-density limit we are considering here, and that 77 = 
d—1 lies for all d in the regime < r\ < 2(d— 1) where 
K acts as an intrinsic ultraviolet cutoff. We conclude that 
at high densities the Coulomb gas is a Fermi liquid in all 
physical dimensions d > 1. The quasi-particle residue 
can be obtained by setting r\ — d — 1 in Eq |2ll I f we 
analytically continue to non- integer dimensions H], we 
see that Z a < 1 for d - 1 < (n/kp)^ 1 . It can be shown 
that r(d, d — 1) ~ | + (3(<i — 1) for small d — 1, so that 
Z Q cx (feir/K^expt-^jy] for d -> 1. 

(e) Quasi-particle damping. If |i? Q | < 00 and the in- 
tegrand in Eq.^(j| is a sufficiently smooth function of q, 
the Fourier integral theorem implies ^"(r, r)| < 00 for 
all r and r, and lim r)T _nxj S a (r, r) = 0. The precise 
way in which S a (r,r) vanishes determines the quasi- 
particle damping. For a conventional Fermi liquid with 
damping 7^ oc (£q) 2 the leading term vanishes for large 
u = max {|r|, |v q t|} as u^ 1 . For interactions of the form 
discussed in (c) we find that S a (r,r) ~ ciu~( d ~ 1 ~ 1 '/ 2 ) + 
C2U~ as u — > 00, where ci and C2 are constants. For 
d < 2 + ^ the first term decays slower than the second 
one, and gives rise to anomalously large corrections. In 
particular, the Coulomb interaction satisfies this condi- 
tion for d < 3, which includes the physically accessible 
case d — 2. The vanishing of S a (r,r) for large u im- 
plies that the effective mass is not renormalized. Thus, 
the definition of £k in our original hamiltonian Hq should 
take the effective mass renormalization into account. 

(f) Relevance of transverse hopping. Experiments 
probing the Luttinger liquid regime are performed on 
anisotropic quasi one-dimensional conductors, consist- 
ing of weakly coupled chains embedded in a three- 
dimensional lattice. If the nature of the three- 
dimensional interaction /™ is such that for vanishing in- 
terchain hopping t± the system is a Luttinger liquid, then 
for any finite t± the quasi-particle residue is finite, so that 
strictly speaking these systems are not Luttinger liquids. 
If tj_ is small compared with the characteristic intra-chain 
hopping energy tii , then the quasi-particle residue is pro- 
portional to (t±/t\\)~ 1 , where the non-universal exponent 
7 > can be identified with the anomalous dimension of 
the Luttinger liquid that would exist for t± — ||l6|| . 

Finally, we would like to point out that our approach 
opens the way for a systematic calculation of corrections 
to the non-interacting boson approximation, which exist 
even in d = 1 if the energy dispersion is not linearized || . 



(n) 

Treating the terms oL' for n > 3 perturbatively within 
one-loop approximation, wc find that for interactions of 
the form f^ a = fae^^^ qc , q c -C kp, the non-Gaussian 
corrections to the bosonized hamiltonian are negligible if 
the dimensionless parameter Aq = Qcj+f^\^\ ^ s small 
compared with unity. Note that Aq vanishes if the energy 
dispersion is linearized, because then the density of states 
is approximated by a constant independent of /x 
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